Instability of black hole formation under small pressure perturbations 
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We investigate here the spectrum of gravitational collapse endstates when arbitrarily small perfect fluid pres- 
sures are introduced in the classic black hole formation scenario as described by Oppenheimer, Snyder and Datt 
(OSD) 1 1]. This extends a previous result on tangential pressures 0] to the more physically realistic scenario 
of perfect fluid collapse. The existence of classes of pressure perturbations is shown explicitly, which has the 
property that injecting any smallest pressure changes the final fate of the dynamical collapse from a black hole 
to a naked singularity. It is therefore seen that any smallest neighborhood of the OSD model, in the space of 
initial data, contains collapse evolutions that go to a naked singularity outcome. This gives an intriguing insight 
on the nature of naked singularity formation in gravitational collapse. 
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I. INTRODUCTION 



One of the most important open issues in the theory and astrophysical applications of modem day black hole and gravitation 
physics is that of the Cosmic Censorship Conjecture (CCC) [iSJ. The CCC postulates that any physically realistic gravitational 
processes must not lead to the formation of a singularity which is not covered by an horizon, thus hiding it from external 
observers in the universe. This of course includes the complete gravitational collapse of a massive star which, if the CCC is true, 
■"sj" must terminate generically into a black hole final state only. Nevertheless in recent years, a wide variety of gravitational collapse 
models have been discovered where the dynamical evolution leads to a naked singularity formation rather than a black hole as 
collapse endstate (see e.g. |01). 

Q The consideration of dynamical evolution is a crucial element of the CCC in the sense that in the original formulation of the 

conjecture it is stated that singularities arising from the dynamical collapse from a regular initial data must be covered by an 
5^ event horizon. Many solutions of Einstein field equations are known which present naked singularities (such as, for example, 
bJO the super-spinning Kerr solutions), nevertheless almost none of these solutions can be obtained as the dynamically evolved final 
state of some initially regular matter configuration. For this reason, over the last decades a great deal of work has been done to 
test the CCC in the few dynamically evolving spacetimes we know. These are typically the scenarios that describe gravitational 
J> collapse in spherical symmetry, and some non-spherical collapse models have also been considered iU. 
\^ In the following, we will consider complete gravitational collapse of a spherical massive matter cloud that leads to the for- 

mation of a strong shell focusing naked singularity at the center In such a case, the super-ultra-dense regions, or the spacetime 
singularity, that forms at the end of collapse would be visible to faraway observers in the universe, rather than being hidden in 
a black hole. Therefore, the question crucial now for the CCC is that of how stable and generic are the naked singularities and 
lO black holes that form in dynamical gravitational collapse of a massive matter cloud. 

' The astrophysical significance of the issue, and the importance of considering the gravitational collapse of a matter cloud 
within the framework of the general relativity theory, with reasonable physical properties for the matter included, stems from 
l the fact that a star more massive than about five to eight times the mass of the Sun, cannot stabilize to a neutron star final state 
J> at the end of its life cycle. It must collapse continually under the force of its own gravity on exhausting its internal nuclear 
fuel, and there are no known forces of nature that would halt such a collapse. General relativity predicts that such a star must 
then terminate into a spacetime singularity where densities and spacetime curvatures blow up and the physical conditions are 
extreme. The CCC assumption that such a singularity is always covered within an event horizon of gravity, is then crucial and 
is at the basis of much of the modern theory and astrophysical applications of black holes today. However, despite the past 
four decades of serious efforts, we do not have as yet available any proof or even any mathematically precise formulation of the 
cosmic censorship hypothesis. 

Actually from the many dynamical gravitational collapse scenarios that have been investigated over the past years, the typical 
conclusion that has emerged is as follows: Depending on the nature of the initial data for the matter cloud in terms of its 
initial density, pressure and velocity profiles from which the collapse evolves, there are dynamical evolutions governed by the 
Einstein equations that take the collapse to either a black hole or naked singularity final state (see e.g. L6J and references therein). 
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Such a naked singularity would be a super-ultra-dense region of extreme gravity that can communicate with faraway observers in 
spacetime, and this hypothetical astrophysical object, if realized in nature, would have radically different observational signatures 
from its black hole counterpart. 

If, however, the occurrence of naked singularities of collapse were special or 'non-generic' in some appropriately well-defined 
sense, then at least the spirit of CCC would be respected. One can impose various sets of physical reasonability and regularity 
conditions under which the collapse of a massive star is to be dynamically evolved to examine its final state. What is really 
needed here is a detailed investigation of the gravitational collapse phenomena within the framework of general relativity, which 
is the only path that can provide useful and adequate insights into the final fate of collapse in terms of either a black hole or naked 
singularity. Investigating how the final black hole state is affected, once a small, general, perfect fluid pressure is introduced 
in the initial data from which the collapse develops would provide a much better understanding of the collapse final states. We 
show here that the models ending in a naked singularity are in fact not 'special', in the sense that they can be generally obtained 
from perfect fluid collapse, where the initial data is arbitrarily close to models leading to a black hole. 

Specifically, we examine here how the evolution of dust collapse models having a black hole final state, is altered when an 
arbitrarily small perfect fluid pressure perturbation is introduced in the matter source. We show explicitly the existence of classes 
of small generic pressure perturbations such that an injection of a small positive (or negative) pressure in the OSD model, or in a 
Tolman-Bondi-Lemaitre (TBL) inhomogeneous dust collapse to a black hole |7], leads the collapse to form a naked singularity, 
rather than a black hole. The classic OSD scenario is the basic paradigm for black hole physics today, and the TBL models 
describe the most general family of dust, i.e. pressureless, collapse solutions. This result is therefore intriguing, because it 
shows that arbitrarily close to the dust black hole model, we have collapse evolutions with non-zero pressures that go to a naked 
singularity final state, thus proving a certain 'instability' of the OSD black hole formation picture against the introduction of 
small pressure perturbations. 

Our method consists in 'injecting' arbitrarily small but generic pressure perturbations in a dynamical dust collapse which was 
originally going to a black hole final state. The pressure is chosen in such a way that it remains small as compared to the energy 
density during the whole collapse and the evolution remains close to the corresponding dust model at all times. We then examine, 
when the small pressures are considered, whether the collapse would evolve to a black hole or a naked singularity. Our analysis 
here shows that in the space of initial data from which the gravitational collapse evolves, any arbitrarily small neighborhood of 
the OSD model would contain collapse evolutions with pressure that go to a naked singularity final fate. While the CCC states 
that the OSD collapse final fate is necessarily replicated for any realistic stellar collapse in nature, the result here shows that an 
arbitrarily small pressure perturbation of the OSD model can change the final outcome of collapse to a naked singularity and 
therefore the OSD black hole may be considered 'unstable' in this sense. Since the pressures within a massive star are very 
important physical forces to take into account, we thus obtain here an important insight into the stability of occurrence of black 
holes as collapse final states. It is such a clarification of the structure of general relativistic collapse, that may provide us a better 
understanding of cosmic censorship, and finally resolve the issue of black hole formation in gravitation theory. 

II. SMALL PRESSURE PERTURBATIONS TO DUST COLLAPSE 

The key feature that characterizes black hole formation in gravitational collapse is, as collapse evolves, trapped surfaces 
and apparent horizon develop at a certain stage within the collapsing cloud, prior to the epoch of formation of the spacetime 
singularity. Then no timelike or null trajectories can escape from the singularity or its vicinity. An event horizon then must form 
hiding the singularity, and the collapse ends in a black hole final state. This is the scenario for an homogeneous pressureless 
collapse model. What we show below is that an arbitrarily small pressure perturbation can radically alter such a scenario. The 
trapped surface formation is then delayed, and this allows the singularity to be naked where families of non-spacelike curves 
can escape from the same. We explicitly identify here such a class of pressure perturbations, but it is by no means the only class 
that can do it. We thus see that the introduction of the slightest pressure in an otherwise pressure free dust model can drastically 
change the final outcome of collapse. 

The most general spherically symmetric metric describing a collapsing matter cloud in comoving coordinates {t, r, 6, (/>) is 
characterized by three free metric functions goo = 9rr = e'^'^'"^*-' and the physical radius of the cloud _R(r, t) = ggg = g^^. 

These are related via the Einstein equations with the energy -momentum tensor, which for perfect fluid matter sources is given 
by Tf — — p; T^' — Tg — Tf — p, where p is matter energy density and p is the pressure in the cloud. The Einstein equations 
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are then. 



F' 



R'^R' 
V' 

P + V 



(3) 



2R' = R'^+R^ , (4) 
F = R{l-G + H) , (5) 

where the dot and prime represent derivatives with respect to t and r respectively and the functions H and G are defined as, 

H = e-2''('^'*)i?2, G = e-2^('-^*)i?'2 . (6) 

Here F is the Misner-Sharp mass of the system, representing the amount of matter enclosed in a radius r at the time t. 

It is known that since the collapse must evolve from a regular initial data, close to the center the radial and tangential pressures 
must be equal for any general collapsing matter field. In fact it can be proven that near the center of the cloud, which is the 
region relevant to our purpose, the pressure gradient must vanish, thus forcing the matter to have a perfect fluid-like behaviour 
18[]. Therefore we work here with a perfect fluid matter cloud as given by the energy-momentum tensor above. 

The model has an additional degree of freedom due to the scale invariance and therefore we can choose the initial time in 
such a way that R{r, ti) = r. We thus introduce the scaling function v{r, t) defined hy R — rv with «(?■, ti) — 1. Collapse is 
described by the condition w < 0, and the singularity is reached at w = 0, where the density diverges. The energy density p is 
regular along the central shell r at any time anteceding the singularity, and that requires the Misner-Sharp mass to have a 
form 

F{r,t)=rH{{r,v{r,t)) , (7) 

with M being finite at the center Also, requiring that the energy density has no cusps at the center and is hence a smooth and 
even function, implies that A/'(0, v) = 0. We then have five equations in six unknowns p, p, M, v, G and v and the system is 
closed if an equation of state relating pressure and density is assumed. In general, however, it is possible and sometimes even 
desirable to study physically valid dynamics satisfying various regularity and energy conditions, without assuming a priori any 
equation of state on which we know little at very high densities of matter 

We therefore choose a suitable physically motivated mass function M as the free function. The physical validity of this choice 
suggests to consider a mass profile arbitrarily close to a well-known collapse scenarios such as the OSD or TBL models. We 
therefore deal with an arbitrarily small pressure perturbation of the TBL model, where by 'small' we mean that the pressure 
remains much smaller than the energy density at all times. 

Once M is fixed, we can evolve the collapse using Einstein equations above and a spacetime singularity develops as collapse 
final state. Such an evolution from regular initial data in terms of initial density and pressure profiles has been studied in detail 
and the conditions for the formation of either a black hole or naked singularity have been worked out |8]. The spacetime 
singularity, corresponding to the epoch t; = 0, is written as the time curve tg (r) which describes the time at which the shell 
labeled by r becomes singular and where v{ts{r), r) = 0. In a neighborhood of the center this is given by, 

t,(r) =to + X2(0)r2+o(r3) . (8) 

The quantity X2 (0) is the tangent to the singularity curve at the origin, and if it has a positive value then the singularity turns 
out to be naked, while otherwise it would be hidden in a black hole. When it vanishes, one has to consider the next order in 
r. Considering collapse of a perfect fluid implies the Misner-Sharp mass F is in general not conserved, so we have to match 
the collapsing cloud with an exterior generalized Vaidya spacetime, which is always possible when the pressure of the matter 
vanishes at the boundary [9]. 

The overall behaviour of the collapsing cloud is determined by the mass function M{r, t), the evolution v{r, t), and the initial 
velocity profiles b{r) for the cloud, which are all not independent and are governed by Einstein equations. The special case of 
homogeneous perfect fluid is obtained when M = M{t), b{r) — k and v — v{t), while the inhomogeneous dust case (TBL) 
is obtained for AI = M (r) (since pr ~ pe, this implies p — 0). Finally the OSD collapse model is obtained when M = Mq, 

h{r) ~ k and v — v{t). 

In the following, for the sake of simplicity, we will consider a constant velocity profile given by b{r) = 1, in analogy with the 
marginally bound collapse in TBL models. 
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We consider an explicit class of perfect fluid collapse models by introducing a small pressure to the TBL scenario, to solve 
Einstein equations in a neighborhood of the center To perturb the TBL model with a small pressure, we must allow in general 

V = v{r,t), rather than v = v{t) only and therefore the simultaneous collapse, which always ends in a black hole, does not 
happen. We take the mass profile to be 

M^Mo + M2iv)r\ (9) 

where Mo is a constant. The pressure perturbation is small when Mq ^| A/2 | at all times. We immediately see that setting 
M2 = C reduces the model to inhomogeneous dust, and A/2 = further gives the Oppenheimer-Snyder-Datt homogeneous dust 
case. To start the collapse from an inhomogeneous dust cloud, triggering the pressure perturbation at a later stage, we take, 

M2iv) = C + e{v) . (10) 

The initial condition A/2(l) = C implies e(l) = 0, with v G [0, 1]. For the present choice of the mass profile, the density and 
pressure, as given by equations ([TJ and ^ are, 

P= jr ,p = PTBL~P+^, — : jrr , (11) 

yz V [v + rv ) 

which are fully determined once the function v is derived from the Einstein equations. 

We now need to determine the behavior of the singularity curve ts{r), and hence to evaluate the quantity X2(0) Ist]. To do 
this, given the choice of A/ above, one can integrate (O in the vicinity of r = 0. Once ly is obtained as a function of M and 

V we can further integrate (3). The function G turns out to be dependent on A/, v and the integration function b{r), which, as 
said before, is related to the velocity of the collapsing shells. Regularity imposes some constraints that are easily satisfied in 
this case. Equation ^ then provides the differential equation for v that constitutes the true equation of motion for the system. 
Integrating it to obtain v{r, t) thus solves the set of Einstein equations. On the other hand, since v is monotonically decreasing 
in time, approaching the singularity u = 0, this equation can be inverted to give the solution as t(r, v), where v is now treated as 
a time coordinate. The singularity curve is then given by <s (r) = t{r, 0), and can be expressed as in equation (O. 

In the present case, taking 6(r) — 1, we can evaluate explicitly X2(0) which turns out to be. 



X.(0) = 4/.-^/., (12) 



with 



h = / {C + Y{v))Z{v)dv , (13) 







h - / W{v)Z{v)dv, (14) 







where we have defined the functions Y, W and Z as 



Yiv) - ' (15) 

W(v) = ev{e + e,^v) , (16) 

A \ -3/2 

4 ev ^ 



We now analyze the condition for the occurrence of naked singularity, namely when X2 (0) is positive, and compare the result 
with the pressureless case. At first we take C = 0, by doing so we are perturbing the OSD model. There are two possibilities 
for the behaviour of the function e, namely e > 0, which implies < and positive pressure, and e < giving e „ > and 
negative pressure. We consider here only the case with positive pressures. Actually, it is easy to see that the negative pressures 
more easily favour the occurrence of naked singularities, therefore it is more useful and physically interesting to check how 
positive pressures affect the black hole formation scenario. 

To see the sufficient conditions that must be required for X2 (0) > 0, we analyze the above two integrals separately. Firstly, by 
evaluating Z. „ we see that if we require. 

Condition la: A/p > max{4t) + e^^vj ,v e [0, 1]} , (18) 
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then Z{v) will be monotonically increasing. We can therefore apply the mean value theorem to Ii, since Y is bounded and 
integrable for v G [0, 1] and Z is bounded, integrable, monotonic, increasing and non negative for v E [0, 1]. Therefore 



h = 



Z{1) I Y{v)dv for some 77 G [0, 1]. 



Now we analyze the behaviour of Y , noting that 1^(0) = e(0) > while 1^(1) = |e_^(l) < 0. Therefore, by Weierstrass 
theorem, we know that there exist v G [0, 1] such that Y{v) = (and Y < Q for v > v, see Fig.l). We shall take e{v) in such a 
way that v is the only zero of Y . Then if 77 > w we have /i < and the first integral in equation (fTZt will be positive. If, on the 
other hand, j] < v then we may always choose e in such a way that 



Condition 2: 



Y{v)dv < 



Y{v)dv , 



(19) 



therefore obtaining again /i < (see Fig.l). 

We turn now our attention to the second integral. We see immediately that it is typically negligible with respect to the first 
one, since it is multiplied by a proportionality factor Trrri- Nevertheless explicitly requiring it to be small implies that we choose 
Mq suitably. We see that to have X2 > we must have 



h 



Again applying the mean value theorem we see that 
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for some cj G [0, 1]. Then X2 will be positive if we choose Afo such that 

Condition lb: > ^ ^^^^ 



(20) 



If we choose e to be at least a cubic function in v of the type e — av^ 



cv + d, it is always possible to choose the four 



parameters a, b, c, d in order to fulfill e(l) = 0, negative and condition 2 (which impUes the values of v and min{y(z;)}). 
Then Mq is chosen in order to satisfy the most stringent between condition la and condition lb. 



black hole 




naked singularity 



FIG. 1: Illustrative plot of e{v) and Y{v) in the case of naked singularity formation with positive pressures. 



We therefore see that it is always possible to find a suitable positive pressure perturbation of the OSD collapse model that, 
however small, will cause the collapse to end in a naked singularity. Furthermore, it can be seen that a perturbation with negative 
pressure can uncover the OSD black hole in the same manner. We emphasize that the above consideration is only to typically 
illustrate what is possible in collapse. The above conditions are sufficient for naked singularity formation but not necessary. In 
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fact, the above example and class given indicates that there might well be other classes of perturbations of the OSD model with 
small perfect fluid pressures, still leading to a naked singularity, and a more detailed discussion will be given elsewhere. 

We note that similar considerations also hold for the pressure perturbations of an inhomogeneous dust TBL model. The 
structure is the same as above, and analogous conclusions hold again, a naked singularity results when we have X2 > 0. For the 
sake of clarity we can see that neglecting the second integral in equation (fT2] l (which we have shown small compared to the first 
one), this is certainly the case whenever C is chosen such that 

C <Di= mm{Y{v),v £ [0, 1]} < . (21) 

On the other hand, values of C such that 

C> D2= max{r(u), u e [0, 1]} > , (22) 

will lead to the formation of a black hole. For C G [Di, D2] the explicit form of e(i>) is what determines the sign of X2 (see 
FigHJ- Again, by putting e = we recover the TBL inhomogeneous dust model, where positive (negative) C leads to black hole 
(naked singularity). It is not difficult to see that a similar reasoning, as given above for the OSD model, applies here. In fact, in 
this case from equation ( fT2] i we have 

1 

X2 = X2\c=o-^J CZ{v)dv , (23) 

with X2\c=a given by the case C = studied above and CZ{v)dv positive whenever C > 0, and negative otherwise. 
Therefore once we evaluate X2\c=q from the procedure discussed above we easily see that all those values of C such that 

^ Jq CZ{v)dv < X2\c=Q will lead to a naked singularity. It is immediate to see that whenever X2\c=a is positive this will 
include positive values of C that for the corresponding TBL model were leading to a black hole. For certain choices of e we 
therefore have models where positive values of C lead to the formation of a naked singularity, even when the TBL collapse went 
to a black hole. Vice versa, for other choices of e we can have models where negative values of C lead to the formation of a 
black hole, whereas the corresponding TBL case was leading collapse to a naked singularity. 

We thus see that the structure of Einstein equations describing gravitational collapse is indeed very rich and complex. As we 
see above, both black hole and naked singularity outcomes are possible, in a general and generic manner in the sense described 
above, evolving from a regular initial data. Of course, in order to be able to check if the black hole that occurs as the endstate of 
the OSD dust collapse is a generic result, we need an unambiguous definition of genericity, which we do not have today. 

In other words, if we call Q the set of physically valid initial data for collapse of a perfect fluid, this will be divided into the 
two possible outcomes as Q — Q^^ |J Q^^ and every point I E Q can be characterized by / = {M{r),p{r), b{r)} (so that 
the OSD initial configuration is given by Iqsd = {Afo, 0, k} G Q^^). We have therefore shown that for every neighborhood 
U{Iosd) C Q, however small, there exist (physically valid) pressure profiles with initial data I <E U such that / G Q'^^ , and 
similarly for the TBL model. In this sense we say that the dust collapse model leading to a black hole is not 'stable' under the 
introduction of small pressure perturbations. Considering the mass and pressure profiles to be expandable in a neighborhood of 
the center we can move from the space of functions in which M and p are defined to the space of the relevant parameters to 
determine the outcome of collapse. This can give some further insights in terms of such parameters. 

This adds to the known results on the OSD model, namely that the introduction of inhomogeneities or a suitable tangential 
pressure perturbation can lead to naked singularity formation [2]. In this sense, we have a strong indication that the OSD 
dust model is 'unstable' in that the initial configurations for its endstates lie on the critical surface separating the two possible 
outcomes of collapse discussed above. 

m. CONCLUDING REMARKS 

We have shown that the introduction of a suitable, though general, arbitrarily small, perfect fluid pressure perturbation in 
the well known OSD collapse model can change the final outcome from a black hole to a naked singularity. The pressure, 
containing only quadratic terms, is chosen in such a way as to satisfy all usual physical requirements such as energy conditions 
and regularity of initial data. 

The model presented above lacks a constitutive relation between the pressure and the energy density. Once such an equation of 
state is introduced the system becomes closed and we do not have the freedom to chose the matter profile at will. Nevertheless, 
assuming an equation of state is not likely to change substantially the picture since the above model is arbitrarily close to the 
well known dust models, and it can be expected to be also arbitrarily close to the solution in the presence an equation of state for 
a suitable choice of the constitutive relation. Further investigation into perfect fluid collapse with an equation of state is currently 
in progress. In this case we still expect the initial configurations in the pressureless case to act as a critical surface separating the 
space of configurations leading to naked singularity from the one leading to black hole. 
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The important question would be whether naked singularities can actually occur in the observable universe and if so, if they 
bear signatures in any way different from the black holes. Indeed, if singularities signal a breakdown of classical gravity when 
very high densities are reached in very small volumes, then these models show that we must consider the possibility that the 
region of spacetime dominated by quantum gravity can affect outside observers and interact with the rest of the universe. 

One should then ask what kind of observational signature these objects bear, if any, and whether such phenomena can possibly 
be observed. The relevant point is if any measurable amount of energy can come out of such ultra-strong gravity regions and, 
if so, in what form. We know that non rotating spherically symmetric collapse, settling to a Schwarzschild black hole, cannot 
emit gravitational waves, still particles and photons can in principle escape the ultra-dense region, thus carrying with them some 
of the energy |10]. Nevertheless, a more accurate description of the phenomenon must take into account rotation and some 
constitutive relation for the matter model. At present this is attainable only within the field of numerical relativity and we hope 
that future research will bring some more light on the nature of the final stages of collapse, since we believe that the next frontier 
of black hole physics will lie in further investigations on the same. 
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